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ABSTRACT 

The turnover in the stellar initial mass function (IMF) at low mass suggests the 
presence of two independent mass functions that combine in different ways above and 
below a characteristic mass given by the thermal Jeans mass in the cloud. In the 
random sampling model introduced earlier, the Salpeter IMF at intermediate to high 
mass follows primarily from the hierarchical structure of interstellar clouds, which 
is sampled by various star formation processes and converted into stars at the local 
dynamical rate. This power law part is independent of the details of star formation 
inside each clump and therefore has a universal character. The fiat part of the IMF 
at low mass is proposed here to result from a second, unrelated, physical process 
that determines only the probability distribution function for final star mass inside a 
clump of a given mass, and is independent of both this clump mass and the overall 
cloud structure. Both processes operate for all potentially unstable clumps in a cloud, 
regardless of mass, but only the first shows up above the thermal Jeans mass, and 
only the second shows up below this mass. Analytical and stochastic models of the 
IMF that are based on the uniform application of these two functions for all masses 
reproduce the observations well. 
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1 INTRODUCTION 

Recent theoretical models suggest that the stellar ini- 
tial mass function may result from star formation in 
hierarchically-structured, or multifractal clouds that are 
characterized by a local conversion time of gas into stars 
that is everywhere equal to the dynamical time on the cor- 
responding mass scale (Elmegreen 1997, 1999a, hereafter pa- 
pers I and II) . This model gives the Salpeter power law slope, 
which was suggested to be the best representation of the 
IMF at intermediate to high mass (see reviews of observa- 
tions in Scalo 1986, 1998; papers I, II, and Massey 1998), 
and it gives a break from that slope to a flattened or possi- 
bly decreasing part at lower mass where interstellar clumps 
cannot become self-gravitating during normal cloud evolu- 
tion. This break point may be at the thermal Jeans mass, 
given by the Bonner-Ebert critical mass 
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for temperature T and total pressure P (thermal + turbulent 
+ magnetic) in the star-forming region. 

The thermal Jeans mass is fairly constant from region to 
region in normal galaxies (Paper II), or at least as constant 
as the IMF is observed to be, considering the limited data 



on the low- mass flattened part. Yet the characteristic mass 
Af J may be varied, if needed, to account for a high mass bias 
in starburst regions (Rieke et al. 1980) or the early Universe 
(Larson 1998) if increases more than P^^^. Such a change 
might be expected in high density regions because they have 
more intense and concentrated star formation. 

A general lack of understanding of the details of star 
formation during the final stages has limited the success of 
the model so far to the power law range, and to all of the 
implications of stochastic, rather than parameterized, star 
formation (Papers I and II). The low mass part of the IMF 
had not been observed well anyway, so the original models 
made no effort to fit any data there. 

There is a growing consensus, however, that the low 
mass IMF becomes approximately flat (slope=0) over a sig- 
nificant range in mass when plotted as a histogram of the 
log of the number of stars per unit logarithmic mass inter- 
val versus the log of the mass. In such a plot, the Salpeter 
(1955) slope is —1.35. This fiattening was suspected over two 
decades ago (Miller & Scalo 1979), but the most recent ob- 
servations are much clearer (Comeron, Rieke, & Rieke 1996; 
Macintosh, et al. 1996; Festin 1997; Hillenbrand 1997; Luh- 
man & Rieke 1998; Cool 1998; Reid 1998; Lada, Lada & 
Muench 1998; Hillenbrand & Carpenter 1999). 
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As a result of these new observations, we suspect that 
the physical origin of the low mass IMF can be understood 
in statistical terms to the same extent as the power law part. 
This paper offers one possible explanation for the flattening 
of the IMF at low mass that may be easily tested with mod- 
ern observations of star formation in resolved clumps. This 
explanation is made here after first reviewing the origin of 
the power law part of the IMF at intermediate to high mass. 



2 THE RANDOM SAMPLING MODEL 

Diffuse interstellar clouds and the pre-star formation parts 
of self-gravitating clouds are generally structured in a hier- 
archical fashion with small clumps inside larger clumps over 
a wide range of scales (see reviews in Scalo 1985; Elmegreen 
& Efremov 2000; for a review of cloud fractal structure, 
see Stiitzki 1999). Stellar masses occur in the midst of this 
clump range, neither at the smallest nor the largest ends, 
and there is no distinction in the non-self-gravitating (pre- 
stellar) clump spectrum indicating where the stellar mass 
spectrum will eventually lie. These observations imply that 
the basic environment for star formation is scale-free, so the 
mass scale for stars has to come from specific physical pro- 
cesses. In papers I and II, we proposed that the mass scale 
arises because of the need for self-gravity to dominate the 
most basic of all forces, thermal pressure, and this leads to 
the M > Mj constraint discussed above. In fact, the break 
point in the power law part of the IMF arises at about Mj 
for normal conditions. 

The Salpeter power-law portion of the IMF was then 
proposed to arise as stars form out of gas structures that 
lie at random levels in the hierarchy (where M > Mj). Per- 
haps the physical process that initiates this is turbulence 
compression followed by gravitational collapse of the com- 
pressed slab (Elmegreen 1993; Klessen, Heitsch, & MacLow 
2000), or perhaps it is an external compression acting on a 
certain gas structure. Regardless, the star formation process 
samples the hierarchical structure of the cloud in this model, 
and builds up a stellar mass spectrum over time. 

This process of sampling where and when a particular 
star forms can only be viewed as random at the present 
time, just as the time and place of rain in a terrestrial cloud 
pattern is random and given only as a probability in most 
weather forecasts. The detailed star formation processes are 
not proposed to be random, only the manifestation of them, 
considering that turbulence is too complex for initial and 
boundary conditions to be followed very far over time and 
space with enough certainty to produce a predictable result. 
Thus we discuss here and in the previous two papers only 
the probability of sampling from various hierarchical levels, 
and we use this probability to generate a stochastic model 
that builds up the IMF numerically after many random tri- 
als. From this point of view, the basic form of the IMF is 
fairly easy to understand, even though the final slope has 
not been derived analytically (i.e., it has been obtained only 
by running the computer model) . 

For example, if the sampling process has a uniform 
probability over all hierarchical levels, then the mass spec- 
trum is M"^dM = M"MlogM exactly (e.g.. Fleck 1996; 
Elmegreen & Falgarone 1996). A dynamically realistic cloud 
will not sample in such a uniform way, however. It will be 



biased towards regions of higher density as they evolve more 
quickly. For most dynamical processes that precede star for- 
mation, including self-gravitational contraction, magnetic 
diffusion, and turbulence, the rate of evolution scales with 
the square root of the average local density. This means 
that lower mass clumps are sampled more often than higher 
mass clumps. As a result, the mass spectrum steepens to 
M~^'^^dlog M . The extra —0.15 in the power comes from 
{D — 3)/2D for fractal dimension D (Paper I). 

In addition to this steepening from density weighting, 
there is another steepening effect from mass competition. 
Once a small mass clump turns into an independent, self- 
gravitating region, the larger mass clump that surrounds it 
has less gas available to make another star. If it ever does 
make a star, then the mass of this star will be less than 
it would have been if the first star inside of it had never 
formed. After these two steepening effects, the IMF becomes 
M~^'^^dlog Af, which is the Salpeter function. 



3 THE FLAT PART OF THE IMF 

The process of random, weighted sampling from hierarchical 
clouds changes below the thermal Jeans mass because most 
of the clumps there cannot form stars at all: there is not 
enough self-gravity to give collapse no matter how the clump 
is put together. As a result, the process of star formation 
stops at sufficiently low clump mass. This does not mean 
that stars smaller than Mj cannot form. Each collapsing 
clump turns into one or more stars with some efficiency less 
than unity, so although the average star mass is proportional 
to the clump mass in this model, the actual star mass can 
be considerably less. There should even be a range of stellar 
masses coming from each self-gravitating clump of a given 
mass, depending on how the clump divides itself into stars 
and how much disk and peripheral gas gets thrown back 
without making stars. 

We have now discovered a curious aspect to this ran- 
dom sampling model when the conversion of clump mass 
into star mass is considered explicitly. That is, the proba- 
bility distribution function for this conversion reveals itself 
clearly below the characteristic mass, and is, in fact, identi- 
cal to the observed form of the IMF there. We also find that 
this clump-to-star distribution function can apply equally 
well above and below Mj, in a self-similar way, but that it 
only shows up below this mass. Thus the power law part 
of the IMF is nearly independent of the details of how a 
clump gets converted into stars (as long as the dynamical 
rate is involved) and depends primarily on the cloud struc- 
ture, whereas the flat part of the IMF depends exclusively on 
the details of clump-to-star conversion and is independent 
of the nature of the cloud's structure. 

To describe this process mathematically, we suppose 
that each self-gravitating clump of mass Mc makes a range of 
star masses Ms such that the probability distribution func- 
tion, P(e), of the relative star mass, e = Ms/Mc, is inde- 
pendent of Mc. This is consistent with the self-similarity 
of star formation that is assumed for the rest of the IMF 
model. The distribution function is written for logarithmic 
intervals as P(e)d log e. The basic point of this paper is that 
P(e) must be approximately constant for all clump masses 
to give the flat part of the IMF at low mass. 
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The final mass function for stars, in logaritlimic inter- 
vals, ns{Ms)d\og Ms, can now be written in terms of the 
mass function for self-gravitating, randomly-chosen clumps, 
nc{Mc)d\og Mc, as 



Us {Ms 



P{e)nc {Ms/e)dloge. 



(2) 



The upper limit to the integral, emax, is the largest relative 
mass for a star that is likely to form from a self-gravitating 
clump. It is perhaps slightly less than unity when Ms » 
Mj, although its precise value is not necessary here. We 
denote it by the constant emax,o and take emax = e-max,o 
when Ms > M,jemax,o- When Ms < Mjernax,o, the efficiency 
can be at most Ms/Mj since the smallest clump that can 
form stars is Mj. Thus 



Ms 



MAX {Mj,Ms/ 
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The lower limit to the integral in equation (g), Emm, is 
not known yet from observations. The ratio of emax,o to e,nin 
will turn out to be the mass range for the flat part of the 
IMF, which is denoted by TZ. It may have a value of ~ 10 or 
more according to recent observations (Sect. I), but future 
observations of lower mass brown dwarfs could extend the 
flat part further and make TZ larger. In any case, we take 



emin — ^jnaxfi/tZ, 



(4) 



Emin does uot depend on whether Ms is greater than or less 
than Mjtmaxfi- 

To solve the integral in equation (P), we convert e to 
Ms/Mc and dloge = — dlogA/c = —dMc/Mc for constant 
Ms from the left hand side. Then, for logarithmic intervals 
in Ms, 



ns{Ms) = Pono 



(5) 



where Pq = P(e) = l//dloge = l/log7^ for all Mc, 
and no comes from the constant in the clump function, 
n(Mc) = noM~^; x — 1.35 from the random sampling 
model and is the same as the power law in the Salpeter 
function. The integral limits are Mc,min ~ Ms /emax = 
MAX {Mj,Ms/ ) and Mc,max = Ms/e 

min — 

MsTl/ emaXjO . 

The solution to equation ^ depends on whether 
Mjemax,o is greater than or less than Ms- For AIs > 

Mjemax,0, 
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For Mj /n<Ms< Mj 
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constant, 



(6) 



(7) 



and for Ms < Mjemaxfi/Ti, this latter result is zero. At 
Ms = Mjemaxfi, the expressions in equations and (Q) 
are the same. A graph of ns{Ms) from these equations is 
shown in figure 1 for IZ = 10, emax,o ~ 1, and Mj — 0.3 
M0. 

Figure 2 shows the result better using the numerical 
model with random sampling described in Papers I and II, 
but now with a stellar mass equal to a random fraction e of 
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Figure 1. Solution to the analytical expression for the IMF from 
equations 6 and 7 using TZ = 10. 



the chosen clump mass. This random fraction is equally dis- 
tributed over a logarithmic range, as specified in the above 
discussion, by using the equation e = TZ^~^ for random vari- 
able r that is distributed uniformly over [0, 1]. Four different 
values of TZ are used in the figure to show how the length of 
the fiattened part equals TZ. The TZ = 1 case shows the pure 
clump selection spectrum, as in papers I and II, but with a 
sharper lower mass cutoff than in the other papers, taken 
here from a failure probability Pf = exp (^— [M^/Mj]''''^ . We 
also took emax,o = 1 for simplicity; the exact value is not 
important (it occurs in the expression for e{r)). For all of 
the cases, there are 10 levels of hierarchical structure in the 
cloud model, with an average of 2 subpieces per piece, and 
an actual number of subpieces per piece distributed as a 
Poisson variable over the interval from 1 to 4. 

The results indicate that for AIs > Mjemax,o, the model 
IMF is a power law with the same power as the clump mass 
spectrum obtained previously. The distribution P (e) does 
not affect the IMF above the break point even though it ap- 
plies there. This is because each chosen clump makes a wide 
range of stellar masses and contributes a flat component of 
width TZ to the local IMF, but the sum of the number of 
stars at each stellar mass is still a power law with the same 
power as the clump spectrum, independent of P (e) . 

The IMF becomes flat below Mjemax.o because all of 
the stars there come from clumps with masses near AIj. 
In fact, the decreasing nature of the clump spectrum to- 
ward higher masses makes clumps with masses very close 
to Mj the favored parents for stars with masses less than 
Mjemax,o- Thus the low mass IMF is determined entirely by 
the separate mass spectrum for stars that form inside each 
self-gravitating clump. 

Evidently, there are two mass spectra for star forma- 
tion, one coming from cloud structure and clump selection, 
giving the Salpeter function, and another coming from star 
formation inside each clump, giving the flat part. The latter 
function actually applies everywhere, but it does not visi- 
bly affect the Salpeter distribution for intermediate to high 
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Figure 2. IMF model with random sampling from a liierarcliical 
cloud for four values of tlie mass range TZ, defined to be the ratio 
of the maximum fraction to the minimum fraction of the clump 
mass that goes into a star. The straight line has a slope of 1.3 

masses because stars there have a wide range (= TZ) of clump 
masses for parents. It only appears for star masses that form 
from the lowest mass clumps that can make a star. 

The model IMF decreases sharply below the smallest 
star mass that can form in a clump of mass Mj, which is 
Mjemaxfi/'R- The parameters tma^.o and TZ, as well as the 
function P(e), depend on the physics of star formation, un- 
like X in the power law part of the IMF, which depends 
primarily on the physics of prestellar cloud structure. 



4 CONCLUSIONS 

The flat part of the IMF is proposed to result from the dis- 
tribution of the ratio of star mass to clump mass, which 
must also be flat in logarithmic intervals. Such a distribu- 
tion shows up only below the mass of the smallest unstable 
clump mass, so in principle, it might only apply there physi- 
cally. However, it could apply equally well for all clumps in a 
star-forming cloud because it does not actually reveal itself 
at masses above this threshold. The low mass flattening and 
eventual turnover below the brown dwarf range depend in 



unspecified ways on the detailed physics of star formation, 
which also contribute to the single characteristic mass, Mj. 
The slope of the power law part of the IMF, above Mj, de- 
pends primarily on cloud structure, although the same com- 
plexities of star formation probably apply there too. This is 
why the Salpeter IMF appears in so many diverse environ- 
ments: the universal character of turbulent cloud structure 
determines the power law nature and the power law itself 
for intermediate to high mass stars, independent of the de- 
tails of star formation. Variations in the IMF from diverse 
physical conditions should be much more pronounced at low 
mass. 

This model may be checked observationally by deter- 
mining the range and distribution function for final star 
mass in resolved clumps of a given mass. The resolution 
of the clumps is important so one can be sure the clumps 
belonging to each individual star or binary pair are mea- 
sured, and not confused with larger clumps that may in- 
clude substructure associationed with separate stars. If this 
model is correct, then the star mass distribution for clumps 
of a given mass will be much fiatter than the Salpeter func- 
tion, and this may be the case for all clump masses, even 
those containing intermediate to high mass stars where the 
Salpeter function applies to the ensemble of stars coming 
from clumps of different masses. 
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